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NONDISPERSIVE DECAY FOR THE CUBIC WAVE 

EQUATION 

ROLAND DONNINGER AND ANIL ZENGINOGLU 



f"') Abstract. We consider the hyperboloidal initial value problem for the 

CN| cubic focusing wave equation 

(-dt +A x )v(t,x)+v(t,x) 3 = 0, i£l 3 . 

Oh 

^M, Without symmetry assumptions, we prove the existence of a co-dimension 

4 Lipschitz manifold of initial data that lead to global solutions in for- 
l/~) ward time which do not scatter to free waves. More precisely, for any 

,- ^ S 6 (0, 1), we construct solutions with the asymptotic behavior 

\\y~V \\LHt,2t)LHB (1 _ S)t ) <*~ 5 + 

^ as t -¥ oo where v (t, x) = ^ an( } _B (1 _^ )t — ji £ I s : \x\ < (1 - <5)£}. 

-I— > 

1. Introduction 
We consider the cubic focusing wave equation 

(-df + A x )v(t,x)+v(t,xf = (1.1) 

.yv in three spatial dimensions. Eq. ( |1.1[ ) admits the conserved energy 

^ E(v(t, -),v t (t, •)) = §||(t;(t, -),v t (t, •))ll^ lxL2(K 3 ) " \\\v(t, -)\\i Hm , 

and it is well-known that solutions with small H 1 x L 2 (R 3 )-norm exist glob- 

r^T ally and scatter to zero [351 ESI GEH E3], whereas solutions with negative 

,__i energy blow up in finite time [19], [27] . There exists an explicit blowup so- 

j>! lution VT(t,x) = ijt^j., which describes a stable blowup regime [11] and the 

blowup speed (but not the profile) of any blowup solution [28] . see also [2] 

for numerical work. By the time translation and reflection symmetries of 



Eq. (1.1) we obtain from vt the explicit solution vo(t,x) = %-, which is 



now global for t > 1 and decays in a nondispersive manner. However, in the 
context of the standard Cauchy problem, where one prescribes data at t = to 
for some to an d considers the evolution for t > to , the role of vo for the study 
of global solutions is unclear because vq has infinite energy. In the present 
paper we argue that this is not a defect of the solution vq but rather a prob- 
lem of the usual viewpoint concerning the Cauchy problem. Consequently, 
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we study a different type of initial value problem for Eq. (1.1) where we 
prescribe data on a spacelike hyperboloid. In this formulation there exists 
a different "energy" which is finite for vq. 

Hyperboloidal initial value formulations have many advantages over the 
standard Cauchy problem and are well-known in numerical and mathemat- 
ical relativity [144 [THj [T71 137] . However, in the mathematical literature on 
wave equations in flat spacetime, hyperboloidal initial value formulations 
are less common (with notable exceptions such as [3]). We provide a thor- 
ough discussion of hyperboloidal methods in Section [2| where we argue that 
the hyperboloidal initial value problem is natural for hyperbolic equations 
in view of the underlying Minkowski geometry. 

To state our main result, we consider a foliation of the future of the 
forward null cone emanating from the origin by spacelike hyperboloids 



E T := (t,i)elxR 3 :i = 4 + yj ^ + M 2 

where T 6 (— oo, 0). Each St is parametrized by 

^:%|C1 3 ^1 4 , <M*)=(- T2 _ T | XrT 2^ X | 2 ) , 

where B R := {X € M 3 : \X\ < R} for R > 0. The ball B ]T \ shrinks 
in time as T — > 0— , but its image under <3?t is an unbounded spacelike 
hypersurface in Minkowski space. The transformation (T, X) i— > <J?t(X) has 
also been used by Christodoulou to study semilinear wave equations in [3] 
and is known as the Kelvin inversion [36]. Note that in four- dimensional 
notation it can be written as X^ i— > ~X ft /(X v X„) (up to a sign in the zero 
component). To illustrate the resulting initial value problem, we plot the 
spacelike hyperboloids St for various values of T £ (— oo, 0) in a spacetime 
diagram (left panel) and in a Penrose diagram (right panel) in Fig. |1.1| 
along with a null surface emanating from the origin. In our formulation of 
the initial value problem we prescribe data on the hypersurface S_i and 
consider the future development. We refer the reader to Section [2] for a 
discussion on hyperboloidal foliations and their relation to wave equations. 
We define a differential operator V„ by 

(V n v) o frrpQ (vo<f> T )(X) 

T 2 — \X\ 2 T 2 — \X\ 2 ' 

which one should think of as the normal derivative to the surface St (al- 
though this is not quite correct due to the additional factor T2 _\ yp )- Ex- 
plicitly, we have 

V n v(t, x) = (t 2 + \x\ 2 )d t v{t, x) + 2tx j djv(t, x) + 2tv(t, x). 
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Figure 1.1. The spacelike hyperboloids £t in a spacetime 
diagram (left panel) and a Penrose diagram (right panel) to- 
gether with the null surface emanating from the origin (thick 



line with 45 degrees to the horizontal). Compare Fig. 2.2 



On each leaf S^ we define the norms 



Il 2 (e t ) 



I if 1 (St) 



Br 



v o $ T (X) 



T 2 -\X\ 



dX. 



and we write 
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$ T (X) = \/2 



t) o $ r (JT) 

X T 2 -|X| 2 

+ |T|" 2 

T 2 - |X| 2 



L 2 (2 T ) 



dX, (1.2) 

We emphasize that 



-T) 



and thus, ||vo||#i(s T ) + llVn^olli^^j,) ~ \T\ 2. Finally, for any subset A C 
M we denote its future domain of dependence by D + (A). With this notation 
at hand, we state our main result. 

Theorem 1.1. There exists a co-dimension 4 Lipschitz manifold Ai of func- 
tions in r7 1 (S_i) x L 2 (£_i) with (0,0) £ Ai such that the following holds. 
For data (/, g) € Ai the hyperboloidal initial value problem 

(-d? + A x )v(t, x) + v(t, x) 3 = 
u|s_! = t>o|s_i +/ 
Vn.li|s_i = V^ols.i +g 

has a unique solution v defined on L> + (S_i) such that 

l T l 5 {h - ^o||ifl(E T ) + ||V«v - V n W || L 2 (ST) ) < |T|3~ 
for all T G [—1, 0). As a consequence, for any 5 G (0, 1), we have 

\\V ~ V \\ L A {t>2t)L A (B{1 _ s)t) <t^ + 
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as t — >■ oo, i.e., v converges to vo in a localized Strichartz sense. 

Some remarks are in order. 

• As usual, by a "solution" we mean a function which solves the equa- 
tion in an appropriate weak sense, not necessarily in the sense of 
classical derivatives. 

• The manifold A4 can be represented as a graph of a Lipschitz func- 
tion. More precisely, let H := U (£_i) x L 2 (S_i) and denote by 
Br(0) the open ball of radius R > around in %. We prove that 
there exists a decomposition ~H = %\ © Hi with dim H2 = 4 and 
a function Fi^n ,6,5(0) ->■ 7i 2 such that M = {u + F(u) : u G 
%in£><5(0)} provided S > is chosen sufficiently small. Furthermore, 
F satisfies 

\\F(u) -F(v)\\ n <S^\\u-v\\ n 

for all u,ve%if\ B 8 (0) and F(0) = 0. 

• The reason for the co-dimension 4 instability of the attractor vq is 



the invariance of Eq. (1.1) under time translations and Lorentz trans- 
forms (combined with the Kelvin inversion). The Lorentz boosts do 
not destroy the nondispersive character of the solution vq whereas 
the time translation does, see the beginning of Section [4] below for a 
more detailed discussion. In this sense, one may say that there exists 
a co-dimension one manifold of data that lead to nondispersive solu- 
tions. However, if one fixes Vq, as we have done in our formulation, 
there are 4 unstable directions. 

There was tremendous recent progress in the understanding of universal 
properties of global solutions to nonlinear wave equations, in particular in 
the energy critical case, see e.g. [131 021 [22] ■ A guiding principle for all 
these studies is the soliton resolution conjecture, i.e., the idea that global 
solutions to nonlinear dispersive equations decouple into solitons plus ra- 
diation as time tends to infinity. It is known that, in such a strict sense, 
soliton resolution does not hold in most cases. One possible obstacle is the 
existence of global solutions which do not scatter. Recently, the first author 
and Krieger constructed nonscattering solutions for the energy critical fo- 
cusing wave equation [9], see also [32] for similar results in the context of 
the nonlinear Schrodinger equation. These solutions are obtained by con- 
sidering a rescaled ground state soliton, the existence of which is typical for 
critical dispersive equations. The cubic wave equation under consideration 
is energy subcritical and does not admit solitons. Consequently, our result 
is of a completely different nature. Instead of considering moving solitons, 
we obtain the nonscattering solutions by perturbing the self-similar solution 
vo(t,x) = ^-. This can only be done in the framework of a hyperboloidal 
initial value formulation because the standard energy for the self-similar 
solution vq is infinite. 
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Another novel feature of our result is a precise description of the data 
which lead to solutions that converge to vq : They lie on a Lipschitz manifold 
of co-dimension 4. In this respect we believe that our result is also interesting 
from the perspective of infinite-dimensional dynamical systems theory for 
wave equations, which is currently a very active field, see e.g. [Ml [2H ES] ■ 

Finally, we mention that the present work is motivated by numerical 
investigations undertaken by Bizoh and the second author [3J. In particular, 
the conformal symmetry for the cubic wave equation has been used in |3j to 
translate the (linear) stability analysis for blowup to asymptotic results for 



decay. We exploit this idea in a similar way: If v solves Eq. (1.1) then n, 
defined by 

u(T,X)- l -f T X ^-"W 



T 2 — \X\ 2 V T 2 — \X\ 2 ' T 2 — \X\ 2 ) T 2 — \X\ 2 

solves (— fly, + Ax)u(T, X) + u(T, X) 3 = 0. The point is that the coordinate 
transformation (t, x) h-> (T, X) with 

t x 



t 2 — \x\ 2 ' £ 2 — |x| 2 

maps the forward lightcone {(£, x) : \x\ < t, t > 0} to the backward lightcone 
{(T, X) : |X| < —T,T < 0} and t — > oo translates into T — > 0— (see 
Fig. |1.1[ ). Moreover, 

T 2 -\X\ 2V ° (~T 2 -|X| 2 ' T 2 -\X\ 2 ) = (-T) =: u o(^,X) 

and thus, we are led to the study of the stability of the self-similar blowup 
solution no in the backward lightcone of the origin. In the context of ra- 
dial symmetry, this problem was recently addressed by the first author and 
Schbrkhuber [TT], see also [TJ El [10] for similar results in the context of wave 
maps, Yang-Mills equations, and supercritical wave equations. However, in 
the present paper we do not assume any symmetry of the data and hence, 
we develop a stability theory similar to [TT] but beyond the radial context. 
Furthermore, the instabilities of no have a different interpretation in the 
current setting and lead to the co-dimension 4 condition in Theorem |1.1| 
whereas the blowup studied in [TT] is stable. The conformal symmetry, al- 
though convenient, does not seem crucial for our argument. It appears that 
one can employ similar techniques to study nondispersive solutions for semi- 
linear wave equations (— df + A x )v(t,x) + v(t, x)\v (t, x)| p_1 = with more 
general p > 3. 

1.1. Notations. The arguments for functions defined on Minkowski space 
are numbered by 0, 1, 2, 3 and we write fl M , u G {0, 1, 2, 3}, for the respective 
derivatives. Our sign convention for the Minkowski metric r/ is (— , +, +, +). 
We use the notation d y for the derivative with respect to the variable y. 
We employ Einstein's summation convention throughout with latin indices 
running from 1 to 3 and greek indices running from to 3, unless otherwise 
stated. We denote by Mj" the set of positive real numbers including 0. 
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The letter C (possibly with indices to indicate dependencies) denotes a 
generic positive constant which may have a different value at each occur- 
rence. The symbol a < b means a < Cb and we abbreviate a < b < a by 
a ~ b. We write f(x) ~ g(x) for x — > a if lint^a l¥l = 1. 

For a closed linear operator L on a Banach space we denote its domain 
by £>(L), its spectrum by <r(L), and its point spectrum by a p (L). We write 
Rl(z) := (z — L) _1 for z G p(L) = C\<r(L). The space of bounded operators 
on a Banach space X is denoted by B(X). 



2. Wave equations and geometry 

In this section, we present the motivation for using hyperboloidal coor- 
dinates in our analysis and provide some background. We discuss the main 
arguments and tools in a pedagogical manner to emphasize the relation be- 
tween spacetime geometry and wave equations for readers not familiar with 
relativistic terminology. 

2.1. Geometric Preliminaries. A spacetime (M,g) is a four-dimensional 
paracompact Hausdorff manifold M with a time-oriented Lorentzian metric 



g. The cubic wave equation (1.1) is posed on the Minkowski spacetime 
(R , rf). In standard time t and Cartesian coordinates (x,y,z) the Minkowski 
metric reads 

rj = -dt 2 + dx 2 + dy 2 + dz 2 , (t, x, y, z) G M 4 . 

Minkowski spacetime is spherically symmetric, i.e., the group SO(3) acts 
non-trivially by isometry on (R 4 ,??). We introduce the quotient space Q = 
R 4 /S*0(3), and the area radius r : Q — > R such that the group orbits 
of points p G Q have area 47rr 2 (p). The area radius can be written as 
r = \J x 2 + y 2 + z 2 with respect to Cartesian coordinates. The flat metric 
can then be written as r\ = q + r 2 da 2 , where q is a rank 2 Lorentzian metric 
and da 2 is the standard metric on S 2 . Choosing the usual angular variables 
for da 2 , we obtain the familiar form of the flat spacetime metric in spherical 
coordinates 

77 = -dt 2 + dr 2 + r 2 (d6 2 + sin 2 6d(f) 2 ) , (t,r, 9, 0) G RxR+ x [0, vr] x [0, 2vr). 

A codimension one submanifold is called a hypersurface and a foliation is 
a one-parameter family of non- intersecting spacelike hypersurfaces. A folia- 
tion can also be defined by a time function from Ai to the real line R, whose 
level sets are the hypersurfaces of the foliation. 

We can restrict our discussion of the interaction between hyperbolic equa- 
tions and spacetime geometry to spherical symmetry without loss of gener- 
ality because the radial direction is sufficient for exploiting the Lorentzian 
structure. Working in the two-dimensional quotient spacetime (Q,q) also 
allows us to illustrate the geometric definitions in two-dimensional plots. 
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2.2. Compactification and Penrose diagrams. It is useful to introduce 
Penrose diagrams to depict global features of time foliations in spherically 
symmetric spacetimes. Penrose presented the construction of the diagrams 
in his study of the asymptotic behavior of gravitational fields in 1963 |34| . 
A beautiful exposition of Penrose diagrams has been given in [Bj. As we 
are working in Minkowski spacetime only, the main features of Penrose di- 
agrams of interest to us is the compactification and the preservation of the 
causal structure. See, for example, [U [21] for the application of Penrose 
compactification to study wave equations in flat spacetime 

The image of the Penrose diagram is a 2-dimensional Minkowski space- 
time with a bounded global null coordinate system. Causal concepts extend 
through the boundary of the map. Consider the rank-2 Minkowski metric q 
on the quotient manifold Q 

q = -dt 2 + dr 2 , (t,r)£txlj. (2.1) 

To map this metric to a global, bounded, null coordinate system, define 
u = t — r and v = t + r for v > u, and compactify by U = arctan u and 
V = arctan v. The quotient metric becomes 

q = / 2TT dVdU, (-tt/2 < U < V < vr/2). 

cos z V cos" 1 U 

Points at infinity with respect to the original coordinates have finite values 

with respect to the compactifying coordinates. The singular behavior of the 

metric in compactifying coordinates at the boundary can be compensated 

by a conformal rescaling with the conformal factor = cos V cos U, so that 

the rescaled metric 

q = Q 2 q=-dUdV, 

is well defined on the domain (—tt/2 <U <V < vr/2) including points that 
are at infinity with respect to q. We say that q can be conformally extended 
beyond infinity. 

The Penrose diagram is then drawn using time and space coordinates 



T = (V + f/)/2 and R = (V - U)/2 (see Fig. |2^2[ ). The resulting metric 
q = —dT 2 + dR 2 is flat. The combined Penrose map is given by 

1 / (T + R\ (T-R 

1 1— > - tan — - — + tan 



2 V V 2 / V 2 

1 / (T + R\ (T-R 

T" i — >■ — { tan — tan 

The boundary dQ = {T = ±(n — R), R G [0, tt]} corresponds to points at 
infinity with respect to the original Minkowski metric. Asymptotic behavior 
of fields on Q can be studied using local differential geometry near this 
boundary where the conformal factor 0, = cosT + cosR vanishes. The part 
of the boundary without the points at R = 0, tt is denoted by J? = {T = 
±(7r — R),R G (0,7r)}. This part is referred to as null infinity because 
null geodesies reach it for an infinite value of their afflne parameter. The 
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Figure 2.2. The level sets of the standard time t depicted 
in a spacetime diagram (left panel) and a Penrose diagram 
(right panel) together with a characteristic line from the ori- 
gin. The boundary of the Penrose diagram includes the spa- 
tial origin and various notions of infinity. Past and future 
timelike infinity are depicted by points i~ and i + . The ver- 
tical line connecting i~ and i + is the spatial origin r = 0. 
Spatial infinity is denoted by the point i°. Null curves reach 
past and future null infinity, denoted by J^~ and J^ + , for 
infinite values of their affine parameter. 



differential of the conformal factor is non-vanishing at J 1 ', dQ\jr 7^ 0, and J" 
consists of two parts J?~ and J?^ referred to as past and future null infinity. 

2.3. Hyperboloidal coordinates and wave equations. Equipped with 
the tools presented in the previous sections we now turn to the interplay 
between wave equations and spacetime geometry. Consider the free wave 
equation 



u tt -Au = -rf'dy.duu = . 



(2.2) 



Radial solutions for the rescaled field v := ru obey the 2-dimensional free 
wave equation 

vtt -v rr = 0, (2.3) 



on (i, r) G 



Ni 



" 



with vanishing boundary condition at the origin. Initial 



data are specified on the t = hypersurface. The general solution to this 
system is such that the data propagate to infinity and leave nothing behind 
due to the validity of Huygens' principle. Intuitively, this behavior seems to 
contradict two well-known properties of the free wave equation: conservation 
of energy and time-reversibility. 

The conserved energy for the free wave equation (12.31) reads 



E(v) 



1 



(v t (t,r) 2 + v r (t,r) 2 )dr. 
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The conservation of energy is counterintuitive because the waves propagate 
to infinity leaving nothing behind. One would expect a natural energy norm 
to decrease rapidly to zero with a non-positive energy flux at infinity. The 
conservation of energy, however, implies that at very late times the solution 
is in some sense similar to the initial state |36| . 

Another counterintuitive property of the free wave equation is its time- 
reversibility, meaning that if u(t,r) solves the equation, so does u(—t,r). 
Data on a Cauchy hypersurface determine the solution at all future and past 
times in contrast to parabolic (dissipative) equations which are solvable only 
forward in time due to loss of energy to the future. 

Both of these counterintuitive properties depend on our description of 
the problem. We can choose coordinates in which energy conservation and 
time-reversibility is violated. Of course, it is always possible to find coor- 
dinates which break symmetries or hide features of an equation. We argue 
below that the hyperboloidal coordinates we employ emphasize the intuitive 
properties of the equation rather than blur them. 

The reason behind the conservation of energy integrated along level sets 



oft can be seen in the Penrose diagram Fig. 2.2 The outgoing characteristic 
line along which the wave propagates to infinity intersects all leaves of the 
t-foliation. When the energy expression is integrated globally, the energy 
of the initial wave will therefore still contribute to the result. The hyper- 
boloidal T-foliation depicted in Fig. |1.1[ however, allows for outgoing null 
rays to leave the leaves of the foliation. Therefore one would expect that 
the energy flux through infinity is negative when integrated along the leaves 
of the hyperboloidal foliation. 



The wave equation (2.3) has the same form in hyperboloidal coordinates 
that we employ 

W T T ~ W RR = , 

where w(T,R) = v(— T ^_ Ri , Ti ^ R2 ). The energy conservation and time- 
reversibility seems valid for this equation as well, but here we have the 
shrinking, bounded spatial domain R £ [0, — T) where T —> 0— . The energy 
integrated along the leaves of this domain 

E{w) = f \ {w T (T, Rf + w R (T, Rf) dR , 
Jo z 

decays in time. The energy flux reads 

f)P 1 

°— = -- (w T (T, -T) - w R (T, -T)f < . 

The energy flux through infinity vanishes only if the solution is constant or 
is propagating along future null infinity. When the solution has an outgoing 
component through future null infinity, the energy decays in time. This 
behavior is in accordance with physical intuition. 

Consider the time reversibility. The equation in the new coordinates is 
time-reversible, but the hyperboloidal initial value problem is not. Formally, 
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this is again a consequence of the time dependence of the spatial domain 
given by R < —T. Geometrically, we see in Fig. 2.3 that the union of the 
past and future domain of dependence of the hyperboloidal surface T = — 1 
covers only a portion of Minkowski spacetime whereas for the Cauchy surface 
t = such a union gives the global spacetime. 





Figure 2.3. Comparison of the future (light gray) and past 
(dark gray) domains of dependence for the Cauchy surface 
t = (left) and the hyperboloidal surface T = — 1 (right). 

In summary, the hyperboloidal foliation given by the Kelvin inversion 
captures quantitatively the propagation of energy to infinity and leads to 
a time-irreversible wave propagation problem. Further, the transformation 
translates asymptotic analysis for t — > oo to local analysis for T — > 0— . 

3. Derivation of the equations and preliminaries 

3.1. First-order formulation and similarity coordinates. We start 
from (— d\ + Ax)u(T,X) + u(T, X) 3 = in the hyperboloidal coordinates 
T = — rs— tzt7> X = a X \„\->. for the rescaled unknown 



t 2 -|x| 2 ' ^ — t 2 -\x\ 2 

u(T,X) 



1 



T 2 -\X\ 



X 



T 2 -\X\ 2, T 2 -\X\ 



As discussed in the introduction, the domain we are interested in is T £ 
[—1,0) and \X\ < \T\. Our intention is to study the stability of the self- 



similar solution uq(T) 
coordinates 



VI 



I- 



Tpy. Thus, it is natural to introduce the similarity 



r = -log(-T), £ 



(-T) 



(3.1) 



with domain r > and |£| < 1. The derivatives transform according to 



d T = e T (d T + ed i >), d 



XJ 



z T d : 



V' 
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This implies 

d\ = e 2T {d 2 T + d T + 2^d^d T + £ J £%% + 2^%) 
and dxjd X j = e 2r d^.dfj. Consequently, for the function 

U(r,0:=n(-e- T ,e~ T 
we obtain from (— d\ + dxjd X i)u(T, X) + u(T, X) 3 = the equation 

[d 2 + 8 T + 2^d v d T - (P k - ^ fc )%% + 2^'%]C/(r,0 = e~ 2T U(T,0 3 . 

In order to get rid of the time-dependent prefactor on the right-hand side, 
we rescale and set [/"(r, £) = e T T/>(r, £) which yields 

[d 2 T + 3d T + 2ed^d T -(^ k -ee)d^d e +Aed^+2]^r,0 = ^{r^f- (3-2) 
The fundamental self-similar solution is given by 

Mr,0 ■= e- T u (-e~ T ,e- T O = V2. 
Writing tp = \/2 + <j) we find the equation 

[d 2 + 3d T + Wdtfdr ~ (6 jk - eC k )d^d e + 4?d$ + 2}<t>(r, 

= 60(r, £) + 3V2<P(t, 2 + (j>(r, 3 - (3-3) 

In summary, we have applied the coordinate transformation 

x 



with inverse 



x 



i-iei 2 ' i-iei 



and 4>{t,£) solves Eq. (3.3) for r > and |£| < 1 if and only if 

v(t,x) = f + W-log ' ,f) (3.4) 



t 2 -|x| 2 ' t 

solves (-d 2 + A x )v(t,x) + v(t,x) 3 = for (t,x) G L> + (£_i). 

We have 8tu(T, X) = e 2r (d T + ^dg + 1)i/j(t, £) and thus, it is natural to 
use the variables <j)\ = (ft, 4>2 = do4> + S, 3 dj(j) + <fi in a first-order formulation. 
We obtain 



i = -£ 3 dj(f)i - 4>i + fa 

dofo = dj&cf)! - gdjfa - 202 + 60i + 3V20 2 + 4>\ . (3.5) 



For later reference we also note that Eq. (3.4) implies 

t 2 d t v(t, x) = -V2- j^ (i-djfa + 0i) + 5^02 (3.6) 

where it is understood, of course, that </>i(t, £) and (J>2(t,£) are evaluated at 
T = -\og w ^ w and£=f. 
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3.2. Norms. Since our approach is perturbative in nature, the function 



space in which we study Eq. (3.5) should be determined by the free version 
of Eq. (3.5), i.e., 

do4>i = -£ J dj4>i -0i + 02 

do4>2 = <9jd J 0i - £ J 3j02 - 202- 

The natural choice for a norm is derived from the standard energy H 1 x I? 
of the free wave equation. In the present formulation this translates into 

II ( M' 7 v)IIhi(b) + II02(t,-)IIl 2 (b) 

where B = {£ G M 3 : |£| < 1}. However, there is a slight technical problem 
since this is only a seminorm (the point is that we are working on the 
bounded domain B). In order to go around this difficulty, let us for the 
moment return to the radial context and consider the free wave equation in 

M 1+3 , 

Utt - Urr - ^Ur = 0, 

in the standard coordinates t and r = \x\. Now we make the following 
observation. The conserved energy is given by 

E(u) = \ \ [u 2 + u 2 r ]r 2 dr. 
Jo 

On the other hand, by setting v = ru, we obtain 

Vtt ~ V rr = 

with conserved energy | ^[v 2 + v 2 .]dr, or, in terms of u, 

/•oo 

E'(u) = \ \ [r 2 u 2 + (ru r + u) 2 ]dr. 
Jo 

The obvious question now is: how are E and E' related? An integra- 
tion by parts shows that E and E' are equivalent, up to a boundary term 
lim r _> 00 ru(r) 2 which may be ignored by assuming some decay at spatial in- 
finity However, if we consider the local energy contained in a ball of radius 
R, the boundary term can no longer be ignored and one has the identity 

i-R r-R 

e 'r(u) ■= \ I [r 2 u 2 + (ru r + u) 2 ]dr = \Ru{R) 2 + \ \ [u 2 + u 2 ]r 2 dr. 
Jo Jo 

The expression on the right-hand side is the standard energy with the term 
2Ru(R) 2 added. This small modification has important consequences be- 
cause unlike the standard energy, this now defines a norm. Furthermore, 
E' R (u) is bounded along the wave flow since it is the local version of a posi- 
tive definite conserved quantity. 

In the nonradial context the above discussion suggests to take 

WM T , ')IIhI(B) + H<M T ' -)\\L 2 (dB) + ||02(T, -)IIl2(B)- 

This norm is not very handy but fortunately, we have equivalence to H 1 x 
L 2 (B) as the following result shows. 
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Lemma 3.1. We haveu 

\m(B) - II/IIh1(B) + Wf\\tf(dB)- 

X 

w 



Proof. For x E M 3 we write r = \x\ and ui = pr. With this notation we have 



/(x) = /(to) and 



i» ( B) = / / I/HI^(«)^ 

'o Jas 



where da denotes the surface measure on the sphere. First, we prove 
WIWl^b) ^ \\f\\m(B) + II / II L 2 (aB)- % density it suffices to consider / E 
C°°(B). The fundamental theorem of calculus and Cauchy-Schwarz imply 

r-r / r-l \ I/ 2 

rf{ruj) = / 9 s [s/(scj)]ds < I / \d r [r f (ru>)}\ 2 dr 



Jo wo 

Expanding the square and integration by parts yield 

\d r [r f (ru)]\ 2 dr = / |5 r /(rw)| 2 r 2 dr + / rd r \f(roj)\ 2 dr 
JO Jo 

+ / |/(TO)| 2 dr 

JO 

= |/H| 2 + / |d r /(TO)| 2 r 2 dr 

JO 



and thus, 

/•l 

|2„2 



r 2 |/(ro;)| 2 < |/(u;)| 2 + / |u;%f(ru;)|Vr/r. 

Jo 



Integrating this inequality over the ball B yields the desired estimate. In 
order to finish the proof, it suffices to show that ||/|| 1,2(35) < ||/||#i(.B), but 
this is just the trace theorem (see e.g. |16j . p. 258, Theorem 1). □ 

4. Linear perturbation theory 

The goal of this section is to develop a functional analytic framework for 
studying the Cauchy problem for the linearized equation 

do4>l = -tfdjfa - 01 + 4> 2 

d <p2 = dj&fa - gdjfa - 20 2 + 60i. (4.1) 

The main difficulty lies with the fact that the differential operators involved 
are not self-adjoint. It is thus natural to apply semigroup theory for studying 



Eq. (4.1 ). Before doing so, however, we commence with a heuristic discussion 
on instabilities. The equation (-9y+Aj)?j(T, X)+u(T, X) 3 = is invariant 

As usual, ||/||l2( SB ) has to be understood in the trace sense. 
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under time translations T t— > T — a and the three Lorentz boosts for each 
direction X 3 

T \— > T cosh a — X 3 sinh a 

X 3 i—T- —T sinh a + X 3 cosh a 

X k ^ X k (k / j) 

where a £ R is a parameter (the rapidity in case of the Lorentz boost). In 
general, if u a is a one-parameter family of solutions to a nonlinear equation 
F(u a ) = 0, one obtains (at least formally) 

= d a F{u a ) = DF(u a )d a u a 

and thus, d a u a is a solution of the linearization of F(u) = at u = u a . In our 
case we linearize around the solution Uo(T, X) = i__ T \ - The time translation 

symmetry yields the one-parameter family u a {T,X) := ^3^ and we have 

d a u a (T,X)\ a= Q = —if?. Taking into account the above transformations 
that led from u to <pi, <p2, we obtain (after a suitable normalization) the 
functions 

Mr,0 = e T , Mr,0 = 2e T (4-2) 



and a simple calculation shows that (4.2) indeed solve Eq. (4.1). Thus, there 



exists a growing solution of Eq. (4.1). Similarly, for the Lorentz boosts we 
consider 

\/2 

a '"' X 3 sinh a — T cosh a 

and thus, d a u a j (T,X)\ a= o = —i^X 3 . By recalling that t^t = S, 3 , this 
yields the functions 

0i(r,O=e j , Mr,0 = ^ j (4-3) 



and it is straightforward to check that (4.3) indeed solve Eq. (4.1). This 



time the solution (4.3) is not growing in r but it is not decaying either. It is 
important to emphasize that in our context, the time translation symmetry 
leads to a real instability. The reason being that u a (T,X) = -^ yields the 
solution 

V *V> X > i2_ N 2 + F _t_^ t + a(t*-\x\*) 

of the original problem. This solution is part of a two-parameter family 
conjectured to describe generic radial solutions of the focusing cubic wave 
equation [3]. If a 7^ 0, v a (t, x) decays like t~ 2 as t — > 00 for each fixed x £ M 3 . 
This is the generic (dispersive) decay. On the other hand, the Lorentz trans- 
forms lead to apparent instabilities since the function u a j yields the solu- 
tion v a> j(t,x) = -. cosh a ^ J sinh - of the original problem which still displays 
the nondispersive decay. Consequently, we expect a co-dimension one man- 
ifold of initial data that lead to nondispersive decay, as mentioned in the 
introduction. Since we are working with a fixed uq, however, there is a 
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four-dimensional unstable subspace of the linearized operator (to be denned 
below). This observation eventually leads to the co-dimension 4 statement 
in our Theorem |1.1[ Note that other symmetries of the equation such as 
scaling, space translations, and space rotations do not play a role in this 
context as the solution uq is invariant under these. 

4.1. A semigroup formulation for the free evolution. We start the 
rigorous treatment by considering the free wave equation in similarity coor- 
dinates given by the system 



n = -£, J dj<pi - 4>i + <t>2 
dofo = dj&fa - tfdjfa - 2(/> 2 - (4.4) 



From Eq. (4.4) we read off the generator 
L u(0 = 



acting on functions in D(L ) := H 2 (B) n C 2 (B\{0}) x H 1 (B) n C 1 (B\{0}) . 



With this notation we rewrite Eq. ( |4.4| ) as an ODE, 

^$(r)=L $(r). 

The appropriate framework for studying such a problem is provided by semi- 
group theory, i.e., our goal is to find a suitable Hilbert space H such that 
there exists a map So : [0, oo) — > B{%) satisfying 

• S (0) = id w , 

• So(t)So(ct) = Sq(t + a) for all r, a > 0, 

• lim T _ !> o+ So(t)u = u for all u G H, 

• lim T _i.o+ 7[So(r)u — u] = Lou for all u G V(Lo) where Lo is the 
closure of Lo. 

Given such an So, the function <3?(r) = So(r)$(0) solves Jr$(r) = L $(r). 
Motivated by the above discussion we define a sesquilinear form on % := 
H\B) n C^B) x L 2 (B) n C(B) by 

(u|v) := I 9j«i(09M0<^+ / «i(w)m(w)d(7(a;)+ I u 2 (0^(0^- 

JB JdB JB 



Lemma 



3.1 



implies that (-| •) is an inner product on Ti and as usual, we denote 
the induced norm by || • ||. Furthermore, we write % for the completion of 
H with respect to || • ||. We remark that W is equivalent to H 1 ^) x L 2 (B) 



as a Banach space by Lemma 3.1 



Proposition 4.1. The operator Lo : T>(Lq) C H. — > % is closable and 
its closure, denoted by Lo, generates a strongly continuous semigroup So : 
[0, oo) — >• B(H) satisfying ||So(t)|| < e _ 2 r for all r > 0. In particular, we 
have <r(Lo) C {z G C : Re z < —\}- 

The proof of Proposition |4.1| requires the following technical lemma. 



Hi 
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Lemma 4.2. Let f G L 2 {B) and e > be arbitrary. Then there exists a 
function u G H 2 (B) n C 2 (B\{0}) such that g G L 2 {B) n C(£\{0}), defined 



by 



tick 



g(o ■■= -{&* - ee)djdMo + sedjuio + f U (o, 



(4.5) 



satisfies \\f - #|k 2 (B) < e - 



Proo/. Since C°°(B) C L 2 (-B) is dense, we can find a g G C°°(B) such that 
11/ — 9lli 2 (B) < §• We consider the equation 



titk 



15 



(^* - ee)d,d k <o + 5^a,-«(o + fu(t) = m) 



(4.6) 



In order to solve Eq. (4.6) we define p(£) = |£|, w(£) = t|t and note that 

Thus, interpreting p and w as new coordinates, we obtain 
g>dju(£) = pdpu(puj) 
^edjdkuiO = ed^ k d^ U (0] ~ Z j d jU (0 = p 2 d 2 p u{ P Lo) 
as well as 



d J dju(puj) 



d ; + 



d 



p 



-<9 P H 



k U j UJ k 



p* 



d^d^k 



d-l 



-Co" 



J a 



^ 



<*;■? 



u(pto) 



where <i = 3 is the spatial dimension. Consequently, Eq. (4.6) can be written 
as 



-(l-p 2 )d p 2 -^ p + 5p^ + f 



£** 



u(puj) = g(pu) 



(4.7) 



where — Ag2 is the Laplace-Beltrami operator on S 2 . The operator — A 5 2 
is self-adjoint on L 2 (S 2 ) and we have <r(— A52) = <r p (— A52) = {£(£ + 1) : 
£ G No}. The eigenspace to the eigenvalue £(£ + 1) is (2£ + l)-dimensional 
and spanned by the spherical harmonics {Y^ m : m £ Z, —£ < m < £} which 
are obtained by restricting harmonic homogeneous polynomials in M 3 to the 
two-sphere S 2 , see e.g. [lj for an up-to-date account of this classical subject. 
We may expand g according to 



9(pu) = Y^9£,m(p) Y i 



l,m\U) 



l,m 



where YITm ^ s shorthand for Xl£o Em=-f ana ^ f° r an y fixed p G [0,1], 
the sum converges in L 2 (S 2 ), see [lj, p. 66, Theorem 2.34. The expansion 
coefficient gn : m(p) is given by 



9e,m(p) = (9(p-)\ Y e,m)L 2 (s 2 ) '■= / g(pu)Yi m (u})da(u) 



s 2 
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and by dominated convergence it follows that g^ m £ C°°[0, 1]. Furthermore, 
by using the identity Y^ m = [£(£ + 1)]~ (— ^.g%)Y^ m and the self-adjointness 
of — A 5 2 on L 2 (S 2 ), we obtain 

9l,m{p) = 7^+i)(5(p-)|(- A S 2 )^,m) L 2( S 2) 
= £(!+!) ((- A S 2 )3(pO|^,m) L 2 (S 2)- 

Consequently, by iterating this argument we see that the smoothness of 5 
implies the pointwise decay ||^, m ||L°°(o,i) — Cm$~ M for any MeN and all 

£ £N. Now we set 

TV 

9n(0 ■= Y^9e,m(\€\)Yi, m (-^) 

and note that ||gjv(p ■) — <Kp ")IIl 2 (5 2 ) -> as JV -> 00. Furthermore, by 
\\g£,m\\L°° (0,1) < ^~ 2 f° r all ^ € N we infer 

sup ||s/vO-)-5(/>-)IIl 2 (.s 2 ) < 1 
pe(o,i) 

for all N £ N and dominated convergence yields 

\\9N ~ 9\\h(B) = / lb/v(p-)-9(/'-)ll! 2 (5 2 )/' 2(i /'^ 

as N — > 00. Thus, we may choose N so large that \\gjy — g\\z?m) < §■ 

By making the ansatz u(puj) = Y^,£ m u e,m(p)Y£ im (u) we derive from Eq. (4.7) 
the (decoupled) system 

- (1 - p 2 )^ - 20, + 5pd p + f + W)] ^, m (p) = %m (p) (4.8) 



for ■£ E No, I < JV, and —£< m < £. Eq. (4.8) has regular singular points at 
p = and p = 1 with Frobenius indices {£, —£ — 1} and {0, — |}, respectively. 
In fact, solutions to Eq. (4.8) can be given in terms of hypergeometric func- 
tions. In order to see this, define a new variable v^ m by u^ m {p) = p e v£ tTn (p 2 ). 
Then, Eq. (4.8) with g^ m = is equivalent to 

z{\ - z)v'l m (z) + [c - (a + b + l)z]v' l>m (z) - abv^ m (z) = (4.9) 

with a = g(|+^), b = a+ |, c= § +^> an d z = p 2 . We immediately obtain 
the two solutions 

A ( v \— Z? /3+M 5+21 3+21 . .A 
J, / 7 \_ Z7 / 3+21 5+21 3.1 _,\ 

where 2-^1 is the standard hypergeometric function, cf. |31| 126] . For later 
reference we also state a third solution, 4>l,ii given by 

Mz) = {l-z)-KF 1 {^^,\;l- z). (4.10) 
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Note that (f>ot is analytic around z = whereas 4>i^ is analytic around z = 1. 
As a matter of fact, <j)\£ can be represented in terms of elementary functions 
and we have 

]__= [(i _ jy^zY^ -(l + VT 



h,e( z ) 



-£- I 

. t- 2 



(4.11) 

see |31j . This immediately shows that |</>i/(z)| — > oo as z — > 0+ which 
implies that ^o,£ and 4>i,e are linearly independent. Transforming back, 
we obtain the two solutions ipj^(p) = p (j)j^(p 2 ), j = 0, 1, of Eq. (4.8) with 
g^ m = 0. By differentiating the Wronskian W(ip 0j t, ipi/) = ipo,e^'i t e-%^i,e 
and inserting the equation, we infer 

^(w,<M'(p) = (i^ - j) w{^^^){p) 

which implies 

W(ik* <M(p) = - - -3 (4-12) 

/r(l — p z )2 

for some constant q. In order to determine the precise value of q, we first 

Tirirp 1" M 3 T 

V^(p) = 2/ +1 ^(p 2 ) + ip e -%Ap 2 )- 

For the following we recall the differentiation formula |31| 

^2-Fi(a, b, c; z) = f 2 F x (a + l,b + l,c + l;z) (4.13) 

which is a direct consequence of the series representation of the hypergeo- 
metric function. Furthermore, by the formula |31j 

lim [(1 - zT+ h - c 2 F x {aAc;z)} = ^^ ^ ~ ° ] , (4-14) 

valid for Re (o + b — c) > 0, we obtain 

■C ^)r(l ) 

4 ^)r(^ 



3+2£\wl- 

lim [(1 - z)Ho, e (z)\ = *" 2 'z?' = ^ (4.15) 



i , », i — j;— M U _* , i 



1 r(^)r(^ 



as well as 



w 3+2l \w3\ 

Jim [(1 - .)%(,)] = r( ^ ) ; ( ^ ) = lh\ 



where we used the identity T(x)F(x + |) = 7T2 2 1 2a T(2x). This yields 
C| = p 2 (l-/) 2 )^(^,^)(/?) 

= P 2 (i - p 2 ) §/<Mp 2 )[V + Vi> 2 ) + ¥-^iAp 2 )} 

- p 2 (l - p 2 )§P^(p 2 )[2/ + Vo,,(p 2 ) + V-Wp 2 )] 
= -2 lim(l-p 2 )t«// 0/ (p 2 ) 

= -2^+5. 
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By the variation of constants formula, a solution to Eq. (4.8) is given by 



ue,m(p) = - 4>o,e(p) 



iPiAp) 



P W(ip 0r t,ipi,i)(s) 1- s 2 



^oA s ) 



9e,m(s) 



ds 



ds. 



/ W(^,^ e )(s) 1- S 2- ( 416 ) 

We claim that U£ >m £ C 2 (0, 1]. By formally differentiating Eq. (4.16) we 
find 



where Ijj, j = 0,1, denote the respective integrals in Eq. (4.16). This 
implies U£ m G C 2 (0, 1) but u" m {p) has an apparent singularity at p = 1. 

We have the asymptotics ipQ£(p)Ii,e(p) — (^ ~ P) an ^ V'i' e(p)IoAp) — 1 as 
p — )■ 1— . Thus, a necessary condition for lim p _>.i_ u'l m (p) to exist is 

a£, m := lim [(1 - p 2 )^ e (p)h,e(p)] = -9i,m0-)- 

This limit can be computed by de l'Hopital's rule, i.e., we write 

hiiP) 



at-, 



lim 



p^i- [(i - p^'Up)]- 1 



lim 
P->1- 



'i» 



(3), 



-[(1 " P 2 )^(P)]- 2 [(1 " P 2 )^» ~ 2P^(P)] 



We have 



lim [(1-p 2 )- 2 l[/p)] 
Pol- 
and thus, it suffices to show that 



\ lim [p 2 i>i,e(p)ge, m (p)] 



fJt, 



Cl 



P-+1 



C( 



1 (1- P 2 )4 3 ](p)-2^(P) 
- — = lim a 

a ^i-(i-p2)3[(l-p2)0)p 



lim 

p->l- 



(l-p 2 )i4 3 j(p)-2p(l-p 2 )l^(p) 
[(l-^)fyUp)] 2 



(4.17) 



Note that 



ip'oiip) = 4p + <PqAp ) + lower order derivatives 



4>) 



^+3j,(3)/ 2 



J 0,£ 



(p ) + lower order derivatives. 



Consequently, from the definition of 4>o,i and Eqs. (4.13), (4.14) we infer 
lim [(1 - p 2 )hoAp)] = 4 lim [(1 - z)hoAz)] = -3c, 

p— >1— z^l— 



,2nL ,(3) 



U(V, 



lim [(1 - p 2 )^»] = 8 lim [(1 - z)2^(z)j 



-15q 
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which proves (4.17). We have g^ m G C°°[0, 1] and thus, in order to prove 
the claim u^ m G C 2 (0, 1], it suffices to show that pi—)- (1 — p 2 )^ t (p)I\^(p) 
belongs to C 1 (0, 1]. We write the integrand in I\^ as 

ip\,i(s) 9i, m (s) \%nm 

where in the following, 0(1) stands for a suitable function in C°°(0, 1]. Con- 

3 ~ 

sequently, we infer I\/(p) = (1 — p) 2(9(1). We have ipoj = a^ipi^ + a,£ipi,e 
where ipi/(p) '■= p l <j>\,t{p 2 ) , see Eq. ( 4.10[ ), and a,£,ae G C are suitable con- 
stants. This yields 

^» = (l-p)-i0(l) + 0(l) 

and thus, (1 - p 2 )ipQ e (p)Ii : e(p) = 0(1) + (1 - p)*0(l). Consequently, p i->- 
(1 — p 2 )ipQj(p)Ii^(p) belongs to C 1 (0, 1] and by de l'Hopital's rule we infer 
Ue, m S C 2 (0, 1] as claimed. 

Next, we turn to the endpoint p = 0. The integrand of I\p is bounded by 
Cep~ +1 and thus, we obtain 

IMp)I<i, ^g{o,i} 

1^1,2^)1 < |logp| 

|/i,Kp)I<p^ +2 , £en,3<e<N 

for all p G (0, 1]. The integrand of Io t i is bounded by Cep e+2 and this implies 

IV(P)I ~ / +3 for a11 P G [M] an d ^ e No, ^ < iV. Thus, we obtain the 
estimates 

i4 fc i(p)i < i 

i^(p)i<p max{1 - M} 

\4 k l(p)\<p 2 - k \iog P \+p 2 - k 

\4%(.P)\Zp 2 ~ k i ^N, 3<^<iV (4.18) 

for all p G (0, 1] and k G {0, 1, 2}. _ 

Now we define the function u : i?\{0} — > C by 

JV 
uCO-X^WAmd^D^mCf)- (4-19) 

From the bounds (4.18) we obtain M \djdku(£)\ < |£| _1 which implies u G 
H 2 (B) Pi C 2 (i?\{0}) and by construction, u satisfies 

where ||/ - £jv||z 2 (b) < 11/ - sIIl 2 (B) + lis - 5aHIl 2 (B) < £• a 



2 Note that Y (u) = -4=. 

' v y v47T 
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Proof of Proposition \4-l\ First note that Lo is densely defined. Further- 
more, we claim that 

Re(L u|u) < -i||u|| 2 (4.20) 

for all u E P(Lo). We write [Lou]^ for the yl-th component of Lou where 
A G {1,2}. Then we have 

d fe [L u]i(£) = -ed 3 d kUl {0 - 2d k ux(0 + dkMO- 

By noting that 



Re [djdkU\d k u\} = \dj[d k uid k u\\ 

?d j f(o=dy[ef(o]-mo 

we infer 

and the divergence theorem implies 



Re / <9 fc [L u]i(£)<9 fc -ui(£)d£ = -\ \ d k ui(uj)d k ui(uj)da{u) 

JB JdB 



I / 9fc«i(03*«i(0<£ 

JB 



+ Re / d k u 2 {i)d k Ul {i)di. 

JB 
Furthermore, we have 



djd J ' ui(£)u2(£)d£ = / u> J djUi(u>)u2(uj)da(oj) 

B JdB 



d jUl (0&u 2 (0dt; 

B 



and 



Re / edMOMOdi = \ [ \u 2 (u)\ 2 da(oj) - § / |u 2 (£)| 2 ^ 

JB JdB JB 

which yields 

Re / [L u] 2 (£)-u 2 (£)<i£ = Re / ui djUi(oj)u2{oj)da{(jj) - \\\ii2\\ 2 L 2 

JB JdB 

-Re / djui(£)diu2(0 d €- U u ^\h(B)- 
In summary, we infer 

Re(L u|u) = -|||«i|||-i (jB) - |||«2||ia( B ) + / A(u)da(u) 



(9B) 
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with 
A(u) = -h\M")\ 2 ~ \Wi(u)\ 2 - llVmMI 2 - \\u 2 (lo)\ 2 



— Re [uj 3 djUi(oj)ui(u)] + Re [uj 3 ' djU\{oS)U2(oj)\ + Re [u 2 (uj)ui(uj)] 

<-\\ Ui {uj)\ 2 

where we have used the inequality 

Re (ab) + Re (ac) - Re (be) < \{\a\ 2 + \b\ 2 + \c\ 2 ), a, 6, c € C 

which follows from < \a — b — c\ 2 . This proves (4.20). 
The estimate (4.20) implies 

|| [A - (L + |)]u|| 2 = A 2 ||u|| 2 - 2ARe ((L + |)u|u) + ||(L + ^)u|| 2 
>A 2 ||u|| 2 

for all A > and u G T>(Lq). Thus, in view of the Lumer-Phillips theorem 
(|15j. p. 83, Theorem 3.15) it suffices to prove density of the range of A — Lo 
for some A > — ^. Let f 6 ?i and e > be arbitrary. We consider the 
equation (A — Lo)u = f. From the first component we infer u 2 = t^djUi + 
(A + \)u\ — /i and inserting this in the second component we arrive at the 
degenerate elliptic problem 

-( < p'*-e j 'e*)5 J -9j fc «(o+2(A+2)^a 7 -«(e)+(A+i)(A+2Mo = /(o (4.21) 

for u = Ul and /(£) := &djfi(Z) + (A + 2)/ x (0 + f 2 (0- Note that by 
assumption we have / 6 L 2 (B). Setting A = \ we infer from Lemma 4.2 the 



existence of functions u G H 2 (B) n C (B\{0}) and g G L 2 (B) such that 

and \\f -g\\ L z(B) < £■ We set u\ := u, u 2 (£) := £ J dju(£) + §«- /i, ffi := /i, 
and 52 (£) := g(£) - ^/i(0 " |/l(0- Then we have u e P(L ), g^, 

l|f - g|| = II/2 - S2||r»(B) = 11/ - &llz»(B) < e 
and by construction, (^ — Lo)u = g. Since f E % and e > were arbitrary, 
this shows that rg (k — Lo) is dense in % which finishes the proof. □ 

4.2. Well-posedness for the linearized problem. Next, we include the 
potential term and consider the system 

do4>\ = -i ] dj4>i - 4>i + 4> 2 

dofo = did* fa - tfdjfa - 2^2 + 60i. (4.22) 

We define an operator L', acting on H, by 

L ' u «» : = ( el 

Then we may rewrite Eq. ( 4.22[ ) as an ODE 
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for a function $ : [0, oo) — > %. 

Lemma 4.3. The operator L := Lo + L' : P(Lo) C % — > % generates a 
strongly-continuous one-parameter semigroup S : [0, oo) — > B(H) satisfying 
||S(t)|| < e'~2 + " '^ T . Furthermore, for the spectrum of the generator we 
have <t(L)\<t(Lo) = cr p (L). 

Proof. The first assertion is an immediate consequence of the Bounded Per- 
turbation Theorem of semigroup theory, see [15], p. 158, Theorem 1.3. In 
order to prove the claim about the spectrum, we note that the operator 
L' : H — > H is compact by the compact ness of the embedding H l (B) <^->- 



L 2 (B) (Rellich-Kondrachov) and Lemma 3.1 Assume that A G cr(L) and 
A £ cr(Lo). Then we may write A — L = [1 — L'Rl (A)](A — Lo). Observe 
that the operator L'Rl (A) is compact. Furthermore, 1 G <t(L'Rl (A)) since 
otherwise we would have A G p(L), a contradiction to our assumption. By 
the spectral theorem for compact operators we infer 1 G o" p (L / Rl (A)) which 
shows that there exists a nontrivial f G T~L such that [1 — L'RL (A)]f = 0. 
Thus, by setting u := Ri, (A)f, we infer u G X>(Lo), u/0, and (A — L)u = 
which implies A G a p (L). □ 

4.3. Spectral analysis of the generator. In order to improve the rough 



growth bound for S given in Lemma 4.3, we need more information on the 



spectrum of L. Thanks to Lemma 4.3 we are only concerned with point 



spectrum. To begin with, we need the following result concerning T>(Lq). 
Lemma 4.4. Let 5 G (0, 1) and u G V(L ). Then 

u\b 1 _ s €H 2 {B 1 „ 5 )xH x {B 1 _ 5 ) 

where B^ s := {£ G M 3 : |f | < 1 - 6}. 

Proof. Let u G D(Lq). By definition of the closure there exists a sequence 
(u n ) C T>(Lq) such that u n — > u and L(ju n — > Lou in H as n — > oo. We set 
f n := L u„ and note that f„ G H X (B) D C\B) x L 2 (B) n C{B) for all n G N 
by the definition of V(L ). We obtain u 2n (£,) = £ j djU ln (£) + itin(£) + /in(£) 
and 

- ($ jk ~ e^djdkUmiO + ^d jUln (0 + 2ui„(0 = /«(0 (4-23) 



where /„(£) := -£ J 0i/i„(O-2/in(O-.MO, cf- Eq. ^- 21 l By assumption 
we have /„—)•/ in L 2 (B) for some / G L 2 (B). Since 

(S jk -eZ k )VjVk>\v\ 2 -\e\ri\ 2 >l\ri\ 2 
for all £ G Bi_g/2 and all ry G M 3 , we see that the differential operator 



in Eq. (4.23) is uniformly elliptic on B 1 _g/ 2 . Thus, by standard elliptic 



regularity theory (see |16j . p. 309, Theorem 1) we obtain the estimate 

WuinWH^B^s) < Cs (Wuinh^B^s^) + \\fn\\L»(B 1 _e/2)) ( 424 ) 
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and since u n — > u in % implies u\ n — > u\ in L 2 (B), we infer u\\b 1 _ s G 
H 2 (Bi_s). Finally, from u 2n (£,) = £ j djUi n (£,) + ui n (£) + /in(£) we conclude 
u 2 \b 1 _ s G H\B^s). □ 

The next result allows us to obtain information on the spectrum of L by 
studying an ODE. For the following we define the space H^ ad (a, b) by 

\\f\W {ab) --= f\np)?p 2 d P + f\f(p)\ 2 P 2 dp. 

Lemma 4.5. Let A G <r p (L). Then there exists an £ G No and a nonzero 
function u G C°°(0, 1) n #^(0, 1) suc/t tfiai 

-(l-pV(p)-^(p) + 2(A + 2WW 

+ [(A + 1)(A + 2) - 6]«(p) + ^u(p) = (4.25) 

/or a// p G (0,1). 

Proof. Let u G £>(L) = £>(Lo) be an eigenvector associated to the eigenvalue 
A G tfp(L). The spectral equation (A — L)u = implies ^2(£) = £ J, 9j"Ui(£) + 
(A + l)tti(0 and 

_(#* _ gftdjdkmiO + 2(A + 2)^'S,-«i(e) 

+ [(A + l)(A + 2)-6]«i(0 = 0, (4.26) 



cf. Eq. (4.21), but this time the derivatives have to be interpreted in the 

weak s 

H\B 1 _ 



weak sense since a priori we merely have u\ G H 2 (Bi_$) n H 1 ^) and «2 £ 
_,j) n L 2 (B) by Lemma 



4.4 



However, by invoking elliptic regularity 
theory ([IB], P- 316, Theorem 3) we see that in fact u% G C co (B)nH 1 (B). As 
always, we write p = |£| and w = tIt . We expand ui in spherical harmonics, 
i.e., 



ui(pui) 



(4.27) 



£,m 

with u^ m (p) = {u\(p •)\Yt,m)L 2 {s 2 ) an< i f° r each fixed p G (0,1), the sum 
converges in L 2 (S 2 ). By dominated convergence and U\ G C°°{B) it follows 
that U£ :m G C°°(0, 1). Similarly, we may expand d p u\{puj) in spherical 
harmonics. The corresponding expansion coefficients are given by 

{d p u l {p-)\Y^ m ) L 2 {S 2 ) = 5 p («i(p-)l^,m)x2(52) = d P ue,m(p) 

where we used dominated convergence and the smoothness of u\ to pull out 
the derivative d p of the inner product. In other words, we may interchange 
the operator d p with the sum in Eq. (4.27). Analogously, we may expand 



A S 2U\(p-) and the corresponding expansion coefficients are 

(A S 2U 1 (p-)\Y £)m ) L 2 {S 2 ) = (u 1 (p-)\A S 2Yt m ) L 2 {S 2 ) = -£(£+l)ut, m (fi). 



Thus, the operator A52 commutes with the sum in (4.27). All differential 



operators that appear in Eq. (4.26) are composed of d p and A52 and it is 



therefore a consequence of Eq. (4.26) that each U£ m satisfies Eq. (4.25) for 



/' 

Jo 
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all p G (0,1). Since at least one u^ m is nonzero, we obtain the desired 
function u G C°°(0, 1). To complete the proof, it remains to show that 
tt/, m € 5^(0,1). We have 

\U£,m(p)\ = \(. u l{p-)\ Y e,m)L2(S 2 )\ < \\ u l(.P IIl 2 (5 2 ) 

and thus, 

l ,i 

\u e , m {p)\ 2 p 2 dp< / ||wi(/0-)lli2( S '2 )J o 2 <i/0 = ||«i||| 2( m- 
o Jo 

Similarly, by dominated convergence, 

\d p Ul, m {p)\ = |(5pUi(p-)l^,m)L2 (5 2)| < ||V«l(p-) 11^2(52) 

and thus, 

l4m(p)lVdP<l|Vui|||2 (B) . 

Consequently, u\ G H^^B) implies u^ m G 3^(0, 1). □ 

Proposition 4.6. -For t/ie spectrum of L we /iaue 

O-(L) c{z G C : Rez < -\} U {0, 1}. 

Furthermore, {0, 1} C 0p(L) and £/ie (geometric) eigenspace of the eigen- 
value 1 is one- dimensional and spanned by 

u(e;i)=(a 

whereas the (geometric) eigenspace of the eigenvalue is three-dimensional 
and spanned by 

u j (Z;0)=(? j Y j G {1,2,3}. 

Proof. First of all, it is a simple exercise to check that Lu(£, 1) = u(£; 1) 
and Luj(£;0) = for j = 1,2,3. Since obviously u(-; 1), Uj(-) 0) G "D(Lo), 
this implies {0, 1} C <r p (L). 

In order to prove the first assertion, let A G cr(L) and assume Re A > — \. 
By Proposition 4.1 we have A ^ o"(Lo) and thus, Lemma 4.3 implies A G 
Op(L). From Lemma 4.5 we infer the existence of a nonzero u G C°°(0, 1) n 
#^(0, 1) satisfying Eq. ( |4.25| ) for p G (0, 1). As before, we reduce Eq. ( |4.25[ ) 
to the hypergeometric differential equation by setting u(p) = p e v(p 2 ). This 
yields 

z{\ - z)v"{z) + [c-(a + b+ l)z]v'(z) - abv(z) = (4.28) 
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with a = \{-l + £ + A), b = \{A + £ + A), c = § + £, and z = p 2 . 
fundamental system of Eq. (4.28) is given byn 



<f>ij(z; A) = 2-Fi (a, &, a + b + 1 — c; 1 — z) 

faj.{z; A) = (1 - z) c - a ~ b 2 F 1 {c -a,c-b,c-a-b + l;l-z) 

and thus, there exist constants q(A) and q(A) such that 

viz) = ct(\)<j>ie(z; A) + c^(A)0i #(z; A). 



A 



The function <^i^(z; A) is analytic around z = 1 whereas (fri/iz; A) ~ (1— z 
as z — > 1— provided A / 0. In the case A = we have (fri/iz; A) ~ log(l — z) 
as z — > 1— . Since u G #^(0,1) implies u G -ff 1 ^,!) and we assume 
Re A > — g, it follows that q(A) = 0. Another fundamental system of 



Eq. (4.28) is given by 

<KeO;A) = 2 -Fi(a,&,c;z) 

o ,£(z; A) = z 1 ~ c 2 i ? i(a - c + 1, b - c + 1, 2 - c; z) 

and since <t>o,i{z; \) ~ z a as z — > 0+, we see that the function p \- > 
P e 4>e,o(p 2 ) does not belong to H^ d {0, §). As a consequence, we must have 
v(z) = de(\)(fto : £(z] A) for some suitable de(\) G C. In summary, we conclude 
that the functions <Pq^{-\ A) and <f>i/(-; A) are linearly dependent and in view 
of the connection formula [31] 



a, < w r(i-c)r(a + 6 + i-c) ^ 

<fii plz A = — — — r(br\p(z A) 

Y ' K ' ' T(a + l-c)T(b + l-cy ' v ' ; 

r( c -i)r(a + b + i-c) I 
+ fRFW M *' A) 

this is possible only if a or b is a pole of the T-function. This yields — a G No 
or -b G N and thus, ±(1 - £ - A) G N or -|(4 + £ + A) G N . The latter 
condition is not satisfied for any £ G No and the former one is satisfied only 
if (£,X) = (0,1) or (£,X) = (1,0) which proves <r(L) C {z G Rez < -|} U 
{0, 1}. Furthermore, the above argument and the derivation in the proof of 



Lemma 4.5 also show that the geometric eigenspaces of the eigenvalues 



and 1 are at most three- and one-dimensional, respectively. □ 

Remark 4.7. According to the discussion at the beginning of Section |4j the 
two unstable eigenvalues 1 and emerge from the time translation and 
Lorentz invariance of the wave equation. 



Strictly speaking, this is only true for c — a — b — —A 7^ 0. In the case A = there 
exists a solution 0i^ which behaves like log(l — z) as z — > 1 — . 
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4.4. Spectral projections. In order to force convergence to the attractor, 
we need to "remove" the eigenvalues and 1 from the spectrum of L. This 
is achieved by the spectral projection 

P:= — f(z-L)- l dz (4.29) 

where the contour 7 is given by the curve 7(5) = \ + |e 27 " s , s E [0, 1]. By 
Proposition |4.6| it follows that "f(s) E p(L) for all s E [0, 1] and thus, the 



[46] it f( 
jq. (|4.29 



integral in Eq. ( |4.29 ) is well-defined as a Riemann integral over a continuous 
function (with values in a Banach space, though). Furthermore, the contour 
7 encloses the two unstable eigenvalues and 1. The operator L decomposes 
into two parts 

L M :rgPnP(L)^rgP, L M u = Lu 

L s : kerPnP(L) -^kerP, L s u = Lu 

and for the spectra we have a(L u ) = {0, 1} as well as a(L s ) = <r(L)\{0, 1}. 
We also emphasize the crucial fact that P commutes with the semigroup 
S(t) and thus, the subspaces rgP and kerP of T~L are invariant under the 
linearized flow. We refer to |2Uj and [12] for these standard facts. However, it 
is important to keep in mind that P is not an orthogonal projection since L 
is not self-adjoint. Consequently, the following statement on the dimension 
of rg P is not trivial. 

Lemma 4.8. The algebraic multiplicities of the eigenvalues 0,1 £ o~ p (L) 
equal their geometric multiplicities. In particular, we have dimrgP = 4. 

Proof. We define the two spectral projections Po and Pi by 
P n = ^/ (z-L^dz, ne{0,l} 

where 70 (s) = \e 27Tis and 71(a) = 1 + \e 2lTis for s E [0,1]. Note that 
P = P + Pi and P Pi = P1P0 = 0, see [20]. By definition, the alge- 
braic multiplicity of the eigenvalue n G o>(L) equals dimrgP n . First, we 
exclude the possibility dimrgP n = 00. Suppose this is true. Then n be- 
longs to the essential spectrum of L, i.e., n — L fails to be semi-Fredholm 
(|20j. p. 239, Theorem 5.28). Since the essential spectrum is invariant under 
compact perturbations, see [20] , p. 244, Theorem 5.35, we infer n G ct(Lq) 



which contradicts the spectral statement in Proposition 4.1 Consequently, 
dimrgP„ < 00. We conclude that the operators L( n ) := L| rg p nrT p(L) are 
in fact finite-dimensional and cr(L( n \) = {n}. This implies that n — L/ n ) 
is nilpotent and thus, there exist m n G N such that (n — L( n )) m " = 0. 
We assume m n to be minimal with this property. If m n = 1 we are done. 
Thus, assume m n > 2. We first consider L( )- Since kerL is spanned 
by {uj(-;0) : j = 1,2,3} by Proposition 4.6, it follows that there exists a 
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u S rgPonP(L) and constants c%, 02,03 G C, not all of them zero, such that 

L (0) u(0 = Lu(£) = £ c,u,(£ 0) = ( £?. 
This implies 1*2 (£) = ^djUiiQ + u i{0 + Cj£ 3 an d thus, 



m=— 1 

As before in the proof of Lemma |4.5[ we expand ui as 



u l(0 = E^,m(|CI)^,m(p 



and find 

- (1 - p 2 )u'{ m {p) - ju' ljTn (p) + Apu' lm (p) - 4u hm (p) + jzu hm (p) = c m p. 

(4.30) 
For at least one m S {—1,0, 1} we have c m 7^ and by normalizing u\^ m 
accordingly, we may assume c m = 1. Of course, Eq. (4.30) with c m = 
is nothing but the spectral equation (4.25) with £ = 1 and A = 0. An 
explicit solution is therefore given by ip{p) = p which may of course also be 
easily checked directly. Another solution is ip(p) := ipo,i{p\ 0) = p4>o t i(p 2 ;0) 
where </>i,o( - ;0) is the hypergeometric function from the proof of Lemma 
We have the asymptotic behavior ip(p) ~ p~ 2 as p — > 0+ and |t/>(p)| ^ 



4.6 



log(l — p)\ as p — >• 1— . By the variation of constants formula we infer that 
ni m must be of the form 



fp ^( s ) s 
Ul,m(p) =C*l>(p) + Cip(p) + ${p) I ^T^Y^^ 

k \ r ^ s ) s a 

-^ p) J Pl wlsjT^ ds 



ds 



(4.31) 



for suitable constants c, c G C, po, p\ € [0, 1] and 



d 



P 2 {l-P 2 ) 

where d G M\{0}. Recall th at u\ £ H l {B) implies ui, m G -£^(0, 1) and by 
considering the behavior of (4.31) as p — > 0+, we see that necessarily 

s 



Vfo) 

pi W(s ) 1 _ S 2 



ds 



which leaves us with 



Ul, m (p) = Clp(p) + lf>(p) j 1 _ ,, 



4>(s) 



ds - i>(p) 



il>(s) 



W(s)l 



;ds. 
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Next, we consider the behavior as p — )■ 1—. Since 

rp ip(s) s 
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po 

1 / 1- 



W(s)l 



;ds 



<1 



for all p G (0, 1) and ip ^ H^ d (^, 1), we must have 



lim , 



jS 4 . Thus, we arrive at a con- 



This, however, is impossible since ^rnx j^2 — a* 
tradiction and our initial assumption too > 2 must be wrong. Consequently, 
from Proposition 4.6 we infer dimrgPo = dimkerL = 3 as claimed. By 
exactly the same type of argument one proves that dimrgPi = 1. □ 

4.5. Resolvent estimates. Our next goal is to obtain existence of the 
resolvent R L (A) G B(H) for A G H_i +e := {z G C : Rez > -\ + e} and |A| 
large. 

Lemma 4.9. Fix e > 0. Then there exists a constant C > such that 
Rl(A) exists as a bounded operator on % for all A G H_i +e with |A| > C. 



Proof. From Proposition 4.1 we know that Rl (A) G B{T~L) for all A G H_ i_ 
with the bound 

I|Rl (A)|| < ' 



ReA + 



i • 



see |15j . p. 55, Theorem 1.10. Furthermore, recall the identity Rl(A) 
R-l (A)[1 - L'Rl (A)] _1 . By definition of L' we have 



L / R L() (A)f 





6[R Lo (A)f] ] 



where we use the notation [gj^ for the fc-th component of the vector g. Set 
u = RL (A)f for a given f G %. Then we have u G £>(Lo) and (A — Lo)u = f 
which implies U2(£) = £, 3 djU\{£) + (A + l)-ui(£) — /i(£), or, equivalently, 



[R Lo (A)f]i(0 



1 



A + 1L 
Consequently, we infer 



^■[R Lo (A)f]i(0 + [R Lo (A)f] 2 (0 + h(0 



[RL (A)f] 



< 



1 



1||L2(B) £ p^y 

+ II/iIIl 2 (b) 
fll 



RL (A)f]i||^ 



(B) 



[R Lo (A)f] 2 ||L2( jB ) 



< 

~ IA + 11 
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which yields ||I/Rl (A)|| < \\+i\ f° r an A ^ H_i +e . We conclude that the 
Neumann series 



[1-L'R Lo (A)]- 1 = ^[L'Rl (A)]' 



fc=0 

converges in norm provided |A| is sufficiently large. This yields the desired 
result. □ 

4.6. Estimates for the linearized evolution. Finally, we obtain im- 
proved growth estimates for the semigroup S from Lemma [4.3| which governs 
the linearized evolution. 



Proposition 4.10. Fix e > 0. Then the semigroup S from Lemma 4-3 
satisfies the estimates 

||S(r)(l - P)f || < Ce^ +e)T \\{\ - P)f || 
||S(r)Pf|| < Ce T ||Pf|| 

for all t > and f G ri. 

Proof. The operator L s is the generator of the subspace semigroup S s de- 
fined by S s (t) := S(r)|k er p. We have cr(L s ) C {z S C : Kez < — -} and 
the resolvent Rl s (A) is the restriction of Rl(A) to kerP. Consequently, 



by Lemma 4.9 we infer ||Rl s (A)|| < 1 for all A G H_i, e and thus, the 
Gearhart-Priiss-Greiner Theorem (see [15] . p. 302, Theorem 1.11) yields the 
semigroup decay ||S s (r)|| < e < -~2+ <E ) T . The estimate for S(r)P follows from 
the fact that rg P is spanned by eigenfunctions of L with eigenvalues and 



1 (Proposition 4.6 and Lemma 4.8). □ 



5. Nonlinear perturbation theory 



In this section we consider the full problem Eq. (3.5), 

docpi = -£ j dj4>i - (f>i + 02 

<9 O 02 = dj&fa - ed j( p 2 - 20 2 + 60i + 3^/20? + 0? (5.1) 



with prescribed initial data at r = 0. An operator formulation of Eq. (5.1) 
is obtained by defining the nonlinearity 



N(u):= (3V24 + n?) 



It is an immediate consequence of the Sobolev embedding H^^B) c -> L P (B), 
p € [1, 6], that N : ri — > ri and we have the estimate 

||N(u)-N(v)||<||u-v||(||u|| + ||v||) (5.2) 
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for all u, v G T~L with ||u||,||v|| < 1. The Cauchy problem for Eq. (5.1) is 
formally equivalent to 

$(0) = u (5.3) 

for a strongly differentiable function <1> : [0, oo) — >■ % where u are the pre- 
scribed data. In fact, we shall consider the weak version of Eq. (5.3) which 
reads 

$(r) = S(r)u+/ S(r - cr)N($(cr))dcr. (5.4) 

Jo 
Since the semigroup S is unstable, one cannot expect to obtain a global 
solution of Eq. (5.4) for general data u G %■ However, on the subspace 
kerP, the semigroup S is stable (Proposition 4.10). In order to isolate the 
instability in the nonlinear context, we formally project Eq. (5.4) to the 
unstable subspace rg P which yields 



P$(t) = S(t)Pu+ / S(T-(j)PN($((T))dff. 
Jo 

This suggests to subtract the "bad" term 

poo 

S(r)Pu+ / S(r - ct)PN($(ct))At 
Jo 

from Eq. ( |5.4[ ) in order to force decay. We obtain the equation 

t*T pOO 

$( T ) = S(t)(1-P)u+ / S(r-<T)N($(<r))d<T- / S(r -a)PN(^(a))da. 

Jo Jo 

(5.5) 

First, we solve Eq. (5.5) and then we relate solutions of Eq. (5.5) to solutions 
ofEq. (IBTl). 



5.1. Solution of the modified equation. We solve Eq. (5.5) by a fixed 
point argument. To this end we define 

K u ($)(t):=S(t)(1-P)u+ f S(t - a)N(${a))da 

Jo 

poo 

- / S(r-<r)PN($(ff))dff. 
Jo 

and show that K u defines a contraction mapping on (a closed subset of) the 
Banach space X, given by 

X := |$ G C([0,oo),?0 : supe^- e ) T ||$(r)|| < ooj 



with norm 



\$\\x ■= supe ( 2 e ) T ||$(r) 



T>0 



where e G (0, ^) is arbitrary but fixed. We further write 

AT* :={$G^: ||$||^ < 5} 
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for the closed ball of radius 5 > in X. 

Proposition 5.1. Let 5 > be sufficiently small and suppose u E % with 
||u|| < 5 2 . Then K u maps Xg to Xg and we have the estimate 

\\K U (*)-K U {*)\\ X <C5\\*-*\\ X 

for all $, \l/ 6 Xg. 

Proof. First observe that K u : Xg ->• C([Q,oo),H) since ||N($(r))|| < 
e (-l+2e)r f Qr any (J, g x s ^ e ^ ave 

/oo 
S(r - o-)PN($(cr))d<7 (5.6) 

which yields 

e^[|«(a)-*((7)[|(||*(a)|| + ||*(a)[|)cto 

/oo 
e (-2+2e)^ a 

< (5e (-l+2£)T|| $ _^|| ;t 

for all $, VlJ G Xg by Proposition |4. 10] On the stable subspace we have 
(l-P)K u ($)(r) = S(r)(l-P)u+ / S(r-a)(l-P)N($(a))da 

JO 

and thus, 
||(1 - P)[K u ($)(r) - K u (W)(r)]|| < £ eH+*)(— )||<&(<t) - *(a)\\ 

X(||$(a)|| + ||*(a)||)da 
< ||$ - $|U<5e ( ~5 +e)r / eK+^a 

<<5eH+^||$-tf||;r 
again by Proposition 4.10| We conclude that 



o 



||K U ($)-K U (*)|U<<5||$-*|U 

for all <3>, $ G A'. By a slight modification of the above argument one 
similarly proves ||K u (<1>)||,y < 5 for all $ £ A^ (here ||u|| < <5 2 is used). □ 



Now we can conclude the existence of a solution to Eq. (5.5) by invoking 
the contraction mapping principle. 

Lemma 5.2. Let 5 > be sufficiently small. Then, for any u £ % with 
|| U H < <5 2 ; £/iere exists a unique solution <3? u G A5 to i±^. (5.5). 



Proof. By Proposition |5.1| we may choose <5 > so small that 
||K U ($)-K u (tf)||*<i||$-*||* 
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for all $, \I/ G X$ and thus, the contraction mapping principle implies the 
existence of a unique $ u G A^ with <1> U = K u (<l> u ). By the definition of K u , 
$ u is a solution to Eq. (5.5). □ 



5.2. Solution of Eq. (5.4). Recall that rgP is spanned by eigenfunctions 



of L with eigenvalues and 1 , see Lemma |4.8| As in the proof of Lemma 
4.8 we write P = Pq + Pi where P n , n G {0, 1}, projects to the geometric 



eigenspace of L associated to the eigenvalue n G a p (L). Consequently, we 
infer S(r)P n = e" r P n . This shows that the "bad" term we subtracted from 
Eq. (5.4) may be written as 



S(r)Pu+ / S(t-ct)PN($ u (<t))At = S(t)[Pu-F(u)] 
Jo 

where F is given by 

l-OO POO 



F(u) :- 







N($ u (<7))d<7-P 



Q 



e- a N($ u (a))da. 



n 



According to Lemma 5.2, the function F is well-defined on B$2 := {u G 
Ti : ||u|| < 5 2 } with values in rgP and this shows that we have effectively 
modified the initial data by adding an element of the 4-dimensional subspace 
rgP of %. Note, however, that the modification depends on the solution 



itself. Consequently, if the initial data for Eq. (5.4) are of the form u + F(u) 
for u G kerP, Eq. (5.4) and (5.5) are equivalent and Lemma 5.2 yields the 



desired solution of Eq. (5.4). We also remark that F(0) = 0. The following 
result implies that the graph 

{u + F(u) : u G kerP, ||u|| < 5 2 } C kerP® rgP 

defines a Lipschitz manifold of co-dimension 4. 



H 



Lemma 5.3. Let 5 > be sufficiently small. Then the function F : Bp 
rg P C 7i satisfies 

||F(u) -F(V)|| < C5||u-v||. 



Xs C X is Lipschitz-continuous. 

Kv(*v)||* 



Proof. First, we claim that u i— > <1> U : Bp 
Indeed, since <3? u = K U (<I> U ) we infer 

||$u - ^v|U < l|K u ($ u ) - K u ($ v )||^ + ||K U ($ V 
< 5\\$ u -$ v \\x + ||u-v|| 
by Proposition |5.1| and the fact that 

||K u ($ v )(r) - K v ($ v )(r)|| = [|S(r)(l - P)(u - v)|| < e^+^\\n 
The claim now follows from ||N(u) — N(v)|| < ||u — v||(||u|| + ||v||). 
We summarize our results in a theorem. 



□ 



Theorem 5.4. Let 5 > be sufficiently small. There exists a co-dimension 
4 Lipschitz manifold Ai C % with G M such that for any u G A4, 



Eq. (5.4) has a solution <1> G Xg. Moreover, $ is unique in C([0,oo),'H). If 
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in addition, u € T>(L) then <3> £ C 1 ([0,oo),'H) and $ solves Eq. (5.3) with 
$(0) = u. 

Proof. The last statement follows from standard results of semigroup theory 
and uniqueness in C([0, oo),%) is a simple exercise. □ 

5.3. Proof of Theorem |1.1| , Theorem |1.1| is now a consequence of Theo- 
rem 



5.4 Eq. (3.4) implies 



and thus, 



v o $ T (X) - vq o $t(X) 

T 2 -\X\ 2 



\V ~ Vq\\ L 2(x t) 



(-T)' 



i og (-r) * ) 



(-T), 



(5.7) 



log(-T),™ 



|T| 



L2(B ) 



= |r|-3[|0 1 (-iog(-r) J .)IU a( B) 

< ITI~ £ 

rs_/ I I 

Similarly, we obtain 

v o $ T (X) -ii o $ T (X) 1 



which yields 



T 2 -|X| 2 



^djM-logi-T),^) 



\ v - v o\\h^h t ) 



T~ 



V0x -log(-T), 



|T| 



L2(B |T |) 



< in 



For the time derivative we infer 
v o $ T (X) - f o o $ T (X) 



9r- 



T 2 -|X| 2 



^ W + ^J^ + <t>) (- log(-T), ^ 



^ 2 (-log(-T), 



(-T) 



and hence, 



|V n u - V„i;o||l2( Et ) = T~ 



2 -log(-T), 



\T\ 



L 2 (fl| T |) 



< ITI 



Finally, we turn to the Strichartz estimate. First, note that the modulus 
of the determinant of the Jacobian of (T, X) t-t (t, x) is (T 2 — |X| 2 )~ 4 . This 
is easily seen by considering the transformation 

X» 

X 11 h+ y» = 

y X a X° 

which has the same Jacobian determinant (up to a sign) since t = —y° and 
x 3 = y° . We obtain 



duV p 



X a X a 5 u >* - 2X U X» 



{x a x°y 



and hence, 



dvy^d^ 



5„ x 
{X a X 



<r\2 
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which yields | det(<9^)| = (X a X a )~ A = (T 2 - \X\ 2 )~ A . Furthermore, note 
that s £ [t, 2t] and x G B(i_^ t imply 

s 



S: = - 
S < - 



9 119 

s z — \x\ z 



> 



It 



it 2 



t 2 - 

1 
< . 

~ 2i 



t > _c&_ 

t 



X 



Consequently, by Eq. (5.7) and Sobolev embedding we infer 



P- W 0|| L 4 (i)2f)L 4 (B ) 



< 



< 



as claimed. 



B (1-*)|S| 
4 



v o $s(X) -vqo $s(X) 



S 2 -\X\ 2 

i,sr'i!-(-io g (-5) 



dXdS 



' \S\)\\L*(B lsl ) dS 



isr'ui-iogi-s^-m^ds 

: :: I " {sf-^ds ~ t" 2+4e 
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